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We use a diagrammatic approach to study low energy physics of a two 
dimensional electron system where the Fermi level is near van-Hove singularies 
in the energy spectrum. We find that in most regions of the ep-T phase 
diagram the system behaves as a normal Fermi liquid rather than a marginal 
Fermi liquid. Particularly, the imaginary part of the self energy is much 
smaller than the excitation energy, which implies well defined quasiparticle 
excitations, and single particle properties are only weakly affected by the 
presence of the van-Hove singularities. The relevance to high temperature 
superconductivity is also discussed. 

Since the discovery of high temperature superconductors (HTSC's), the role played by 
a saddle point in the energy spectrum, known as a van-Hove singularity (vHS), has been 
actively debated in the literature The experimental observation of vHS's near the 

Fermi surface in many HTSC's || has led to the so called van-Hove scenario to explain the 
normal state properties of these materials 0-0], as well as the superconducting properties 
In this scenario, the presence of the vHS's, as well as weak electron electron (e-e) 
interactions, are used to explain the main body of experimental data. 

One of the more notable properties of the normal state is the linear temperature de- 
pendence of the resistivity |J. This has been explained within the van- Hove scenario by 
arguing that when the temperature is larger than the energy difference between the Fermi 
level and the vHS's, the e-e scattering rate is linear in temperature More generally, 

in this temperature region the imaginary part of the self energy has been claimed to have 
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the behavior ImE (k, e) oc max(e, T) for a large range of energies e, thus explaining many 
other normal state properties ||. From a theoretical view, the linear energy dependence of 
ImE (k, e) differs from the quadratic dependence of a regular Fermi liquid (FL), prompting 
some workers to classify the system as a marginal Fermi liquid (MFL) || . 

Recently, it has been suggested |7J that the situation is more complicated, and that in 
a certain region of the ejr-T phase diagram (e^ is relative to the vHS's), FL theory breaks 
down altogether and a non Fermi liquid (NFL) exists instead. The resulting phase diagram 
is shown in Fig. [I], and can be summarized as follows: When €f is small enough, or T is low 
enough (specific criteria will be given further on), the system behaves as a NFL. In other 
regions of the phase diagram the system behaves as either a regular FL or a MFL: When 
T < \ep\, the affect of the vHS's is not felt and ImE oc T 2 , indicating a regular FL behavior. 
In the opposite case the temperature is high enough so that quasiparticles near the Fermi 
level feel the affect of the vHS's, resulting in the linear temperature dependence of ImE 
characteristic of a MFL. 

In the present work we use a standard diagrammatic approach to study the low energy 
physics of a system of weakly interacting electrons with the Fermi level near vHS's. We 
show that the leading order perturbation calculation of the self energy |3|-|5|] is invalid, and 
instead we need to renormalize the self energy by summing diagrams in the particle-particle 
channel. Dzyaloshinskii |7]] used a more complex version of this scheme to show the existence 
of a NFL in the bottom left corner of the phase diagram in Fig. [l|. However, we concentrate 
rather on the implications for the Fermi liquid properties in the region not occupied by a 
NFL. For this region, our summation scheme is strictly valid, and our main result is that the 
energy dependence of ImE (k, e) is always weaker than linear, even in the part of the phase 
diagram previously thought to contain a MFL. Furthermore, we show that the quasiparticle 
properties, namely the effective mass and the quasiparticle weight, are only weakly affected 
by the presence of the vHS's. 

These results imply that the van-Hove scenario in it's present form cannot explain fully 
the normal state properties of HTSC's. These include the linear temperature dependence 
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of the conductivity ||, the suppression of the quasiparticle weight ||, and the extended 
nature of most vHS's ||. Theoretically, the results are important since ImS (k, e) -C e, so 
that quasiparticles are well defined excitations, as required by FL theory. Thus, we conclude 
that regular FL theory is valid in all parts of the phase diagram not occupied by the NFL. 

The basic system we discuss here is that of a square lattice with an arbitrary energy 
spectrum. Due to symmetry, any band of the energy spectrum will have two equal energy 
saddle points per Brillouin zone (located at wave vectors (0,7r/a) and (vr/a, 0), a being the 
lattice constant). When these points are close to the Fermi level, we shall see that the low 
energy physics of the system is dominated by electrons occupying regions in k-space in their 
vicinity. This enables us to describe the system using the model Hamiltonian 

n — C k C k,<T C k,<r 25 C k+q,a C k'-q,cr' C k',cr' C k,cr 
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Here S is the sample area, and c£l (cj^) creates (annihilates) electrons with spin a and 
wave vector k relative to saddle point i{= 1,2). The spectrum near these points is 

J 1 ) _ f l2 _ f p (2) _ . , 2 _ . ,2 (2 \ 

where t x = h 2 /2m x , t y = h 2 /2m y and m x ,m y are the effective masses. If \t x —t y \ 
t x ,t y , there is said to be significant nesting between the two saddle point. This is not the 
general case however, so throughout this work we consider only the non-nested case, meaning 
\t x — t y \ ~ t x ~ t y . g (g) is a k-independent coupling constant representing processes in 
which the initial and final saddle points of a given electron are the same (different). Strictly, 
each of these constants should be split into two, depending on whether or not both incoming 
electrons are at the same saddle point. However, we will show that only processes where 
both incoming electrons are at the same saddle point are relevant, thus allowing us to ignore 
this point. To complete the definition of our model we assume a wave vector cutoff k c ~ 1/a, 
which defines an energy cutoff e c ~ max (t x , t y )k 2 c . 
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To understand the basic physics of this model we first consider a simpler model with only 
one saddle point instead of two. Thus, we omit the g coupling term in the Hamiltonian, 
and have only one species of electrons whose spectrum is given by e k in Eq. (Q). For this 
model it is convenient to perform a change of variables to e k = tk x k y , where t = y/t x t y . 
Our original assumption regarding t x and t y means that t ~ t x ,t y and that the new form 
of the spectrum also has a cutoff ~ k c . The density of states (DOS) for this spectrum is 
D(e) « (\/2txH) In (tk^/e) rs (l/27r 2 t) In (e c /e), which diverges at small e, leading to the 
so called vHS. This means that the main contribution to the DOS comes from states lying 
close to the saddle point, which justifies our original assumption that the low energy physics 
is dominated by such states. Furthermore, a simple estimate shows that an electron with 
energy e<e c will typically have a wave vector k satisfying 

ln ifL^i n ifL^i n i£ >L ( 3 ) 

I ™x I I I I H 

The above equation, which means that all logarithmic factors appearing in our model are 
large and of the same order, is central to this work and will be used throughout. 

We begin our treatment of the simpler model by considering leading order contributions 
to X (k, e), which come from the diagrams in Fig. |2|a and ^b. When |e| < \ep\, it is easy 
to show that these diagrams give regular FL behavior, meaning ImS oc e 2 . Instead, we 
concentrate on the more interesting case |e| > \ep\. Also, for simplicity we neglect finite 
temperature effects, which serve only to replace e with ~ T when |e| < T. Under these 
conditions, we obtain 
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~ 9 6 ^ln2 + (l-ln2)^W(-) for lei > |e k | 



Im£(k,e) = 

l> 8vr 3 t 2 V'"" ' v ~ e J "* V e 

The main point here is that the on shell self energy is linear in the energy, which is a signature 

of MFL behavior. Using the analytic properties of the self energy, we may calculate it's real 

part. Since this is logarithmically large, it's main contribution is determined by ImE (k, e) 

with lei ^> led, so that 
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Re£(k,e) =_^[(ln2) e + (l-ln2) ek ]ln 2 ( ) . (5) 

8nH z \m&x(\€F\, \e\) / 

We note that this result is equivalent to that obtained by Dzyaloshinskii [0], who calculated 
ReS (k, e) directly. Also, it is a further indication of MFL behavior, since it means that the 
quasiparticle properties are strongly renormalized in the low energy limit. 

We next consider higher order contributions to the self energy. It has been shown JTTJ 
that all such contributions may be accounted for by correcting the bare Greens functions 
of diagrams |2|a and or by replacing one of the interaction lines with a vertex. We will 
show shortly that corrections to the bare Greens functions are small, thus leaving us with the 
vertex replacements. The possible vertices are shown in Fig. ||]c-^]e, and an example of a third 
order diagram including such vertices is shown in diagram [|f. The singular nature of these 
vertices has been discussed many times within the context of the van- Hove scenario [|]||||. 
It has been shown that the particle-particle channel has a square logarithmic divergence, i.e. 
T c ~ t~ l g 2 ln 2 (e c /e), whereas the particle-hole channel only diverges as a single logarithm, 
i.e. ~ T e ~ t^g 2 ln(e c /e) (here e is the scale of the incoming energies of the vertices, and 
the subscripts c, d and e refer to Fig. |2|). Thus, at low enough energies, we may neglect the 
particle-hole channel, and need only consider corrections to the self energy coming from the 
particle-particle channel, T c . Furthermore, for energies such that (g/t) ln 2 (e c /e) > 1, these 
corrections are larger then the leading order contribution to the self energy, thus signaling 
the breakdown of leading order perturbation theory. 

The above arguments can easily be extended to higher order diagrams, which means 



that the self energy must be renormalized by summing the particle-particle ladder JL1]] and 
inserting it into diagrams §a and |2|b, as shown in Fig. |^. For |e| < \ep\ we still obtain 
ImS (k, e) oc e 2 , however, the important effect of the renormalization can be seen in the 
opposite case, when |e| > \ep\- Here we obtain 
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ImE (k, e) 
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l + (g/MH) ln 2 (e c /e) 



(6) 



^ln2 + (1 - In 2)— ) 7~~~i~\ lei » |e k | 
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which is the main result of our work. From it, it is clear that as long as (g/t) ln 2 (e c /e) > 1, 
the on shell self energy satisfies ImE (k, e) ~ e/ln 4 (e c /e) <C e. Thus even when |e| > \ep\ 
Landau quasiparticles are well defined, and the system behaves as a regular FL rather than 
a MFL. Using analyticity, we obtain for the real part of the self energy 

ReE(k, £ ) =^Kln2), + (l-. n 2, fk ]l„(l„(— ^)). (7) 

Comparing to Eq. ([|), we see that the real part of the self energy is strongly renormalized, 
so the affect of the vHS is much weaker then in the leading order case. Using Eq. (^) 
to calculate corrections to the quasiparticle weight and effective mass, we see they are ~ 
(g/t) ln(ln(e c /e)). Since our renormalization procedure is valid for the region (g/t) ln 2 (e c /e) > 
1, these corrections can generally be considered small. Furthermore, in most samples \ep\ is 
finite (see below), so the correction are small even for very low energies. This justifies our 
not including corrections to the bare Green's function in our renormalization scheme. 

However, it is clear that the energy £, defined approximately by (g/t) ln(ln(e c /£)) ~ 1, 
is an energy scale below which our approximation becomes invalid. As long as |ep| > £, 
we may apply our results all the way down to T = 0. When |ep| < £, we may assume 
our results describe the system as long as T > £. This estimate for the breakdown of our 
approximation is consistent with the work of Dzyaloshinskii 0. There, by renormalizing 
the Greens function in a self consistent manner, it was shown that the quasiparticle weight 
becomes identically zero at an energy defined by (g In 2/tn 2 ) ln(ln(e c /£)) = 1. This signifies 
the onset of the NFL in the bottom left hand corner of Fig. [I]. Thus, we conclude that our 
results remain valid as long as FL theory remains stable, i.e. as long as max (T, |e_p|) > £. 

At this point it is instructive to consider typical values of e c /e for the Cuprates. Since 
the energy distance of the vHS's from the Fermi level is typically 10 — 20 mev ||, and e c ~ 1 
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ev, we obtain € c j\ep\ ~ 50 - 100, so the factor ln 4 (e c /e) appearing in the denominator of 
ImS can be as large as ~ 100. Furthermore, due to the high power of the logarithm, the 
dependence of ImS on e should easily be distinguishable from a linear dependence. Thus, 
the linear temperature dependence of the resistivity cannot simply be explained by the 
presence of the vHS's. On the other hand, for the factor (g/t) ln(ln(e c /e)) appearing in ReS 
to be large, we would need say, ln(ln(e c /e)) > 2 (assuming weak to intermediate coupling). 
This in turn means that € c /\e\ > 10 3 , or |e| < 1 mev. Even in samples with very small 
\cf\, such an energy scale would probably not be observable due to c-axis coupling, or other 
smearing mechanisms. Therefore, the presence of the vHS is unlikely to explain the strong 
renormalization of quasiparticle properties observed in the Cuprates. 

Next we consider the full model given by Eq. ([!]). Besides the intra saddle point processes 
already discussed, we now include coupling between the saddle points (inter saddle point 
processes). An example of a leading order inter saddle point contribution to the self energy 
is shown in Fig. ||]g. For the non-nested case we consider here, it can be shown that all 
such contributions are smaller than the intra saddle point ones discussed earlier ||. What 
of higher order inter saddle point contributions, which can be constructed by inserting inter 
saddle point vertices into diagrams ^|a and ^|b? Of the various inter saddle point vertices, only 
the particle-particle vertex shown in Fig. |2|h has the necessary ln 2 (e c /e) energy dependence 
to compete with vertex 0c. Thus, we construct the modified ladder approximation shown in 
Fig. |], and obtain for the renormalized scattering vertex 

r 9 + ff 2 «A 2 + g 2 aA 2 

~ 1 + gaA 2 + gaA 2 + (gaA 2 ) 2 - (gaA 2 ) 2 ' 1 ' 

where A = ln(e c /e), and a = 4/ ((27r) 2 t^. We see that as long as g > g we still have 
T oc A -2 for large A, and all our previous results for the self energy remain valid to within 
prefactors of the logarithmic factors (which now depend on g as well as g). If g < g then the 
denominator of Eq. (Q) becomes zero for large enough A, leading to a BCS type instability 
first pointed out by Gonzalez et si ||. Here we do not dwell on the question of whether 
or not g < g, we merely note that g represents processes with large momentum transfer, 



whereas g represents processes with small momentum transfer. Therefore it is reasonable to 
assume that g > g and that no instability develops. 

In conclusion, we have shown that in all regions of the phase diagram where Fermi liquid 
theory remains valid, the system behave as a regular FL, rather than a MFL. This is based 
on the fact that ImE (k, e) < e even in the region previously thought to be occupied by 
a MFL. Furthermore, the effect of the vHS's on the quasiparticle properties, reflected by 
the real part of the self energy, has been shown to be weak. Our work implies that the 
van- Hove scenario, in it's present form, cannot explain all the normal state properties of 
high temperature superconductors. 

These results differ not only from the leading order perturbation results 0-0, but also 
from some more recent attempts to go beyond leading order perturbation theory [f||12j. The 
renormalization procedure used is strictly valid as long as FL theory remains valid, meaning 
we are not near the non Fermi liquid region of the phase diagram in Fig. |l|. For this region, 
the more complex theory due to Dzyaloshinskii [7[ should be used. 
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FIG. 1. A schematic phase diagram, as reflected by previous works. Here we show that in the 
region labeled as a MFL, the system actually behaves as a regular FL. tp = means that the 
Fermi level coincides with the energy of the vHS's. 

FIG. 2. Diagrams used in this work. Solid (double) lines represent electron propagators near 
the first (second) saddle point, whereas dashed (double) lines represent g (g) interaction lines. 
Diagrams (a)-(f) refer to both the full model and the single saddle point model, whereas diagrams 
(g) and (h) refer only to the full model, (a) and (b) are leading order diagrams for the self 
energy; (c), (d) and (e) are vertices used to contruct higher order self energy diagrams; (f) is an 
example of such a third order diagram; (g) is an example of a leading order inter saddle point self 
energy diagram; and (h) is the only inter saddle point vertex which needs to be included in the 
renormalization of the self energy. 

FIG. 3. The ladder approximation for the self energy in the single saddle point model. V is the 
renormalized scattering vertex. 

FIG. 4. The ladder approximation for the full model. T is the renormalized vertex corresponding 
to g, whereas T corresponds to g. 
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